Giant spin orbit interaction due to rotating magnetic fields in graphene nanoribbons 
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We theoretically study graphene nanoribbons in the presence of spatially varying magnetic fields 
produced e.g. by nanomagnets. We show both analytically and numerically that an exceptionally 
large Rashba spin orbit interaction (SOI) of the order of 10 meV can be produced by the non- 
uniform magnetic field. As a consequence, helical modes exist in armchair nanoribbons that exhibit 
nearly perfect spin polarization and are robust against boundary defects. This paves the way to 
realizing spin filter devices in graphene nanoribbons in the temperature regime of a few Kelvins. 
If a nanoribbon in the helical regime is in proximity contact to an s-wave superconductor, the 
nanoribbon can be tuned into a topological phase sustaining Majorana fermions. 

PACS numbers: 73.22.Pr, 75.70.Tj, 73.63.Fg, 72.25.-b 



I. INTRODUCTION 

The last decade has seen remarkable progress in the 
physics and fabrication of graphene-based systems!^ 
The recent advances in producing graphene nanoribbons 
(GNRs) enable to assemble them with well-defined edges, 
in particular of armchair typeP^ Moreover, it has been 
shown that the presence of adatoms can significantly in- 
crease the strength of the spin orbit interaction (SOI) 
of Rashba typeP All this together makes nanoribbons 
promising candidates for spintronics effects. In partic- 
ular, generation of helical states, modes which trans- 
port opposite spins in opposite directions, is of great 
interest. Such modes were pr opo sed in semiconducting 
nanowiresP carbon nanotubesf^ bilayer graphene P and 
experimentally reported for quantum wires in GaAs hole 
gases.^ The y fi nd applications in spin-filters^ Cooper 
pair splittersp^and, i n contact with an s-wave supercon- 
ductor, they provide a platform for Majorana fermions 
with non-abelian braiding statistics.^ 

In the present work we propose a novel way to gener- 
ate a giant effective SOI in GNRs by spatially varying 
magnetic fields that can be produced by nanomagnets P2I 
This approach has an advantage over using adatoms be- 
cause the surface of graphene is not in tunnel- contact 
with other atoms, which usually leads to high disorder 
with strong intervalley scattering. As we will see, large 
values of SOI result in helical modes of nearly perfect 
polarization. Moreover, nanoribbons, in stark contrast 
to semiconducting nanowires, have considerably larger 
subband splittings, allowing for a superior control of the 
number of propagating modes and of the gaps that are 
characteristic for the helical regime. 

Further, our proposal is a next step in bringing topo- 
logical features to graphene systems. Topological states 
proposed by Kane and Mele^l turned out to be experi- 
mentally undetectable due to the small intrinsic SOI of 
graphene. In contast, we show here that if a GNR in 
the helical regime is brought into proximity to an s-wave 
superconductor, the system can be tuned into a topolog- 
ical phase that supports Majorana fermions. This opens 



up the possibility to use GNR for topological quantum 
computing. 

The low-energy physics of armchair GNRs is character- 
ized by broken valley degeneracy enforced by the bound- 
ary effects.^ To generate helical states we also need to 
lift the spin degeneracy. This can be achieved by mag- 
netic fields in two ways: by a uniform magnetic field and 
Rashba SOI or by a spatially varying magnetic field. The 
chemical potential should be tuned inside the gap opened, 
leading to a helical regime. We will study these two sce- 
narios both analytically and numerically. Moreover, we 
will show numerically that the presence of helical states 
is robust against small non-idealities of the GNR edges. 
This shows that our proposal is realistic and experimen- 
tally feasible. 



II. GRAPHENE NANORIBBON 

GNRs are strips of graphene, a honeycomb lattice de- 
fined by translation vectors ai and &2 and composed of 
two types of atoms, A and B, see Fig. [Tj The GNR axis 




FIG. 1. An armchair GNR formed by a finite strip of graphene 
aligned along the z-axis and of width W in i-direction. The 
GNR is composed of two types of atoms A (blue dot) and 
B (green dot) and is characterized by hexagons in real space 
with translation vectors ai and &2. The low-energy physics 
is determined by the momenta k = (k x ,k z ) around the two 
valleys K = — K' = (47r/3o, 0). Nanomagnets (green slabs) 
placed with period A„ on the sides of the GNR provide a 
spatially varying magnetic field B n (red arrows). 
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is chosen along the z-axis and has a finite width in x di- 
rection. GNRs are usually characterized by a width W 
and a chiral angle 9, the angle between the GNR axis and 
ai. We only consider armchair nanoribbons for which 9 
is equal to n/2. 

Graphene can be analyzed in the framework of the 
tight-binding approach. The effective Hamiltonian in- 
cludes hoppings of electrons between neighboring sites, 

H = Uj,\\>c\ x c jX >- (1) 

<ij>,A,A' 

Here, Ci\ are the standard electron operators, i and 
j are nearest-neighbor sites, and A, A' are spin projec- 
tions on the z-axis. Without SOI, the spin is conserved 
and the hopping amplitude becomes £y,AA' = hj^XX't 
where Uj is spin-independent. It is more convenient to 
treat Hq in momentum space (k x ,k z ). The low-energy 
physics of graphene is determined by two valleys around 
K = — K' = (47r/3a,0), where a = |ai| is the lattice 
constant. Wavefunctions can be represented in the form 
ip = J2t<j ^Te' Tjfl1 , where t = ±1 corresponds to K/K' 
and a = ±1 to the A/B sublattice. The Hamiltonian for 
the slowly- varying wavefuctions <$> ar {x,z) is written in 
terms of the Pauli matrices <Ji (r,-), acting on the sublat- 
tice (valley) degrees of freedom, as 

Hq = hv F (T 3 k x ai + k z a 2 ). (2) 

Here, k z (k x ) is the longitudinal (transverse) momen- 
tum calculated from a Dirac point, and vp is the Fermi 
velocity. From now on we work in the basis $ = 

(4>AK, 4>BK, 4>AK',(i>BK<)- 

A GNR, in contrast to a graphene sheet, is of finite 
width W leading to well-gapped subbands. In order to 
impose open boundary conditions on a GNR consisting 
of N unit cells in transverse direction, we effectively ex- 
tend the GNR by two unit cells, so that the width is 
equal to W — (N + 2)a and impose vanishing bound- 
ary conditions on these virtual sites, tp(0, z) = and 
ip(W , z) = Op^ This leads to quantization of the trans- 
verse momentum k x , K x + k x = 7m/ W , where n is an 
integer. 

If the width of the GNR is such that N = 3M + 1, 
where M is a positive integer, the GNR is metallic.^ The 
low-energy spectrum is linear and given by E a — "fhv F k z , 
where the isospin 7± 1, corresponding to the eigenvalues 
of the Pauli matrix a?,. The higher levels are two-fold 
degenerate (apart from spin, see below) and gapped by 
Trhv F /3(M + l)a (see Fig.^). 

If the GNR width is such that N = 3M (or N = 
3M + 2), where M is a positive integer, the GNR 
is semiconducting.^ The gap at k z — is given by 
2hv F \k x nm \, with |fc™ m | = n/3(N + 2)a . In case 
of a semiconducting GNR all orbital states are non- 
degenerate (see Fig. [5J:). This means that the bound- 
aries induce intervalley coupling and break the valley 
degeneracy^Nill 




FIG. 2. The spectrum of an armchair GNR obtained by 
numerical diagonalization of the tight-binding Hamiltonian 
Ho +H ao +H Z . The low-energy spectrum is linear for metallic 
[(a,b) with N = 82] and quadratic for semiconducting [(c,d) 
with N = 81] GNRs. The SOI [(b) A so = 1 meV and (d) 
A so = 5 meV] lifts the spin degeneracy, so the spectrum con- 
sists of (b) two Dirac cones or (c) two parabolas shifted by 
k so from zero (shown by green lines). For a metalllic GNR 
each branch is characterized not only by the spin projection 
s but also by the isospin 7. The solid (dashed) lines corre- 
spond to 7 = 1 (7 = — 1), see (b). While for a semiconducting 
GNR a magnetic field, Az = 0.1 meV, alone opens a gap 2A 9 
[(d)], we also need to include intervalley scattering (modeled 
by fluctuations in on-site energies) for the metallic GNR [(b)] . 
If the chemical potential /1 is tuned inside the gap, the sys- 
tem is in the helical regime with nearly perfect polarization, 
(s x ) w 0.99, in both cases. 

If the chemical potential fj, crosses only the lowest level 
of the spectrum, there are two states propagating in op- 
posite directions with opposite isospins a. However, so 
far we have not taken spin into account, which will lead to 
four states at the Fermi level in total. As we will see next, 
this degeneracy can also be lifted if we include Rashba 
SOI and a uniform magnetic field or, equivalently, a spa- 
tially varying magnetic field. 



III. RASHBA SPIN ORBIT INTERACTION 

The Rashba SOI arises from breaking inversion sym- 
metry. This can be caused by an electric field E ext ap- 
plied perpendicular to the GNR plane, or alternatively, 
by adatoms, which produce local electric fields. In the 
first case, the S OI is quite small, A so = eE ext £, with 
C = 4x 1(T 5 nirPliaifor realistic fields E ext - 1 V/fim. 
In the second case, the strength of the SOI is sig- 
nificantly increased by doping, and values for A so of 
10 — 100 meV have been observecP. The general form 
of the Rashba SOI Hamiltonian can be obtained from 
symmetry considerations,^ 

H so = A so (t 3 s z <7i - s x a 2 ), (3) 
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where the Pauli matrices Si act on the spin. 



IV. ROTATING MAGNETIC FIELD. 

An alternative approach to generate helical modes is to 
apply a spatially varying magnetic fieldPsHU Such a field 
can be produced by nanomagnets^ or by vortices of a su- 
perconductor. We emphasize that this scheme requires 
not perfect periodicity of the field but just a substan- 
tial weight of the Fourier component at twice the Fermi 
wavevector. Moreover, this mechanism is valid for both, 
rotating and linearly oscillating fields. For simplicity we 
assume in this section that the Rashba SOI discussed 
above is negligible. The case when both, a spatially vary- 
ing magnetic field and Rashba SOI, are present was dis- 
cussed recently in the context of nanowires^ and shown 
to lead to a number of striking effects such as fractionally 
charged fermionsP^l 

First, we consider a field rotating in a plane perpen- 
dicular to the GNR, leading to the Zeeman term 



H„ = &z [s v cos(k n z) + s z sin(fc n z)] 



(4) 



where A z = g/isB/2, fis is the Bohr magneton, g the 
g-factor, and X n = 2i:/k n the period of the rotating field. 
It is convenient to analyze the position-dependent H amil- 
tonian Ho+H„ in the spin-dependent rotating basisl^SlSH 
The unitary transformation U n = exp(ik n zs x /2) brings 
the Hamiltonian back to a GNR with in-plane Rashba 
SOI and uniform field perpendicular to the GNR, 

H 1 - = Ul(H + H^)U n = H a + A z s y + A n so s x a 2 . (5) 

In a similar way, a field rotating in the plane of a GNR, 



H\ = A z [s x cos(k n z) + s z sin(fc„z)] , 



(6) 



is equivalent to a GNR with out-of-plane Rashba SOI 
together with a uniform field along the x-axis, 

HW = Ul(H + Hl)U n = H + A z s x + A n so s y a 2 . (7) 

The induced SOI favors the direction of spin perpendicu- 
lar to the applied rotating magnetic field, and its strength 
is given by A" Q = hvpk n /2, independent of the amplitude 
A Z - For example, A" D is equal to 10 meV for nanomag- 
nets placed with a period of 200 nm. 



V. HELICAL MODES. 

The spectrum of H^- (or by analogy of H") can be 
easily found using perturbation theory. Taking into ac- 
count that realistically Az -C A" c , we treat the Zee- 
man term as a small perturbation. The induced SOI, 
given by A™ s x a 2 , leads to spin-dependent shifts of the 
fc 2 -momenta by k so = A" Q /hvp — k n /2, both for the 
metallic and the semiconducting GNRs, see Figs. [2b and 



[2ji. Every level is characterized by the spin projection 
s = ±1 on the x-axis, so the spin part of the wavefunc- 
tions, | a), is an eigenstate of the Pauli matrix s x . The 
corresponding spectrum and wavefunctions that satisfy 
the vanishing boundary conditions (for ip) are given by 



a*,*. = e <*(*.+-*«)(_i 7) i J f 7> _i)| s ) J ( 8 ) 

Ey,s = jfrv F (k z + sk so ) (9) 



for a metallic GNR and 



- 1 r,k 7 ' 



E±, s = ±hv FV / (k™ in ) 2 + (k z + sk so ) 2 



)\s), (10) 
(11) 



for a semiconducting GNR. Here we use the notation 
e i Va = [fc ™„ _ - (fcz + skso )]/ y J{k^Y + {k z + S k so f. 

A uniform magnetic field that is perpendicular to 
the spin-quantization axis defined by the SOI results 
in the opening of a gap 2A g at k z — 0. Using the 
wavefunctions given by Eq. (10), we can show that 
A, 



Azk™ m ■ 



Az for a semicon- 



ducting GNR. The spin polarization in this state is 
given by \{s x )\ « 1 - (A z fc™ m /4hv F k 2 so ) 2 . In contrast 
to that, a metallic GNR possesses an additional sym- 
metry. Each branch is characterized not only by spin 
(s = ±1) but also by isospin (er = ±1), see Fig. [2Jd. 
Thus, a magnetic field alone cannot lift the degeneracy 
at k z = 0. However, if we include also terms breaking 
the sublattice symmetry, such as intervalley scattering 
described by Hkk 1 = Akk'Ti, a gap will be opened. 
Here, Akk i is the strength of the intervalley scatter- 
ing, which can be caused by impurities or fluctuations 
in the on-site potential. Assuming Az,Akk' "C A^, 
the gap becomes 2A g — 2Akk'Az/A so in leading or- 
der. The spin polarization of the helical states is given 
by \(s x )\ wl- (A z /A so ) 2 . We note that for both semi- 
conducting and metallic GNR, Az limits the size of the 
gap A g . 

We note that H is equivalent to the Hamiltonian de- 
scribing a GNR in the presence of Rashba SOI and a uni- 
form magnetic field applied in perpendicular y-direction, 
H tot = Hq + H so + H z [see Eqs. and ^} in first 
order perturbation theory in the SOI. Here, the Zeeman 
term is given by Hz = Azs y . The wavefunctions given 
by Eqs. (JsJ) and ( 10 1 are eigenstates of the Pauli matrix 
Ti, so the diagonal matrix element of T3 is zero. This 
leads to the result that the term T^,s z a\ in the Rashba 
Hamiltonian H so averages out in first order perturbation 
theory, and H to t is indeed equivalent to H . This means 
that the effect of the SOI is a spin-dependent shift of k z 
by k so = A so /hvp. Similarly, the uniform magnetic field 
opens a gap at k z — 0, which can be as big as 10 K for a 
field of about 10 T. 

An alternative approach to above perturbation the- 
ory is to analyze the GNR with Rashba SOI an- 
alytically. For graphene the spectrum of the ef- 
fective Hamiltonian Hq + H so is given by Ej_± — 
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± (A so + jy/(hv F k x ) 2 + (hv F k z ) 2 + A 2 so j , where the 

index j is equal to 1 (—1) for the highest (lowest) electron 
level, and the ± sign distinguishes between electrons and 
holes. The SOI lifts the spin degeneracy, however, the 
valley degeneracy is maintained, and T3 is a good quan- 
tum number. Analogously to RefP^, we search for a sum 
over the eigenstates ip T , q (x) of H + H so , ip E ' k *(x,z) — 



YlrqbT,q''pT, , q"( x )i sucn tnat tne boundary conditions 
are satisfied. The index q = (j, ±) distinguishes be- 
tween four wavevectors satisfying E,j t ±{k x = ±fcj-) = E, 
hvpki.2 = \JE 2 — (hvFk z ) 2 ± 2_EA so . We also intro- 
duce new variables 6 and 7, via cos 9 = hvpk z /E and 
-\/2 sin 6> sin 7 = hvpki/E. We allow for real as well as 
imaginary values of k\ y 2, 9, and 7. The spectrum of a 
metallic GNR is then given implicitly by 



tan 



k ? W 



k x W 



koW 



sin(27) 



sin (kiW) sin (k2W) sin(27). 



(12) 



The exact solution defined by Eq. ( 12 ) can be analyzed 
analytically by means of Tailor expansion. For example, 
if A so <C Hvpk z , we get E — ±hvpk z ± A so , which is in 
agreement with previous perturbative calculations. 



VI. NUMERICS 




To check our analytical results numerically, we extend 
the tight-binding Hamiltonian H Q by allowing for hop- 
pings with spin-flip, 



<ij>,\,\' 



SAA'CjA' + 



H.c. 



(13) 



in such a way that H so is equivalent to the Rashba SOI in 
the low-energy sector. Here, saa' is a vector composed of 
the Pauli matrices, and spin-dependent hopping elements 
are defined as Uij = — (3A so /4)z x e,j. A unit vector eij 
points along the bond between two sites i and j. The re- 
sults of the numerical diagonalization of the Hamiltonian 
Hq +H so +Hz are presented in Fig. [2] where the Zeeman 
term corresponding to a magnetic field B is modeled as 



H z = ^ <4 A B • saa'Qa' 

i,A,A' 



(14) 



As shown in Fig. |2j the numerical results fully confirm 
the analytical calculations. 



FIG. 3. Defects on the edges of a metallic armchair GNR 
(TV = 82) result in opening of a gap at zero energy. In the nu- 
merical diagonalization, the defects were modeled by omitting 
two atoms on the edges, which was assumed to be periodic 
with period Id = 5\/3a- We see that the spectrum changes 
only little and the qualitative features of a metallic armchair 
GNR are maintained. 



GNRs ! 16 * 17 l We emphasize that for the scenario of heli- 
cal modes developed above we do not need any specific 
symmetries. Thus, our proposal is robust against edge 
defects. 

The scenario with a rotating magnetic field is even 
more universal. The only criterion is that the Fermi 
wavevector kp is not too large, typically kp/K x should 
be smaller than 1CP 2 . This is a natural limit resulting 
from the fact that the period of rotation of a magnetic 
field should be much larger than the lattice constant. 



VIII. MAJORANA FERMIONS 



VII. STABILITY AGAINST EDGE DEFECTS. 

The spectrum of GNRs is known to be sensitive to the 
specific form of the edges. For example, the linear spec- 
trum of a metallic GNR becomes parabolic for non-ideal 
armchair boundaries (see Fig. [3]). In contrast to that, 
subband gaps are only slightly modified for semiconduct- 
ing GNRs. We conclude that the valley degeneracy, in 
general, is lifted due to strong intervalley mixing induced 
by the boundaries and this is a property of all armchair 



Next, we consider GNRs in the helical regime brought 
near to a superconductor. If /1 is tuned inside the gap 
opened by the field, the two propagating modes are he- 
lical. The proximity-induced superconductivity in the 
GNR leads to the coupling between such states and gaps 
the spectrum. The system can be effectively described in 
the linearized model for the exterior (x = e : states with 
momenta close to the Fermi momentum, k e — kp) and 
the interior branches (x — h states with nearly zero mo- 
menta, ki = 0)P^The electron operator is represented as 
^(z) = Yj P =±i x =e i e vfexZ * PX , where the sum runs over 
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the right (R, p = 1) and left (L, p 
effective Hamiltonian becomes 



1) movers. The 



H = -ihv F p 3 xsdx + —rV3pi{l + xs) + A s i] 2 p 2 , (15) 



where the Pauli matrices x% {Vi) ac t m t ne interior- 
exterior branch (electron-hole) space, and W = 



fee^Le, *L. *L> *«. *«i *L *D- Following 
Refsp2122l ; we find that the criterion for the t opological 
phase transition is given by A g > \x 2 + A^. In terms 

Af/A^ 
for a metallic 



of Zeeman energy this gives Az > A so y / /? 
(semiconducting) GNR. 



IX. CONCLUSIONS 

We have shown that helical modes can be generated 
in graphene nanoribbons by a spatially varying magnetic 
field or by Rashba spin orbit interaction with a uniform 
magnetic field. We have demonstrated that the open- 
ing of the gap is universal for both semiconducting and 
metallic graphene armchair nanoribbons independent of 
the mechanism that induces the spin orbit interaction, 
leading to a helical regime with nearly perfect spin po- 
larization. Moreover, we have checked numerically that 
the helical regime is robust against boundary defects. All 
this makes graphene nanoribbons promising candidates 
for spin effects and spintronics applications. 
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